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Abstract 

We studied cross nuclear and mitochondrial gene pools of brown hares (Lepus europaeus) from 
three local populations in Britain and two in New Zealand, to test the hypothesis of reduced genet¬ 
ic variability in hares from New Zealand resulting from few founders originating from Britain. Multi¬ 
locus allozyme electrophoresis of 52 protein loci and analysis of restriction fragment length poly¬ 
morphisms of total mitochondrial DNA based on 16 hexanucleotid-cleaving restriction enzymes 
were employed in 119 and 36 hares, respectively. Observed and expected average heterozygosities, 
rates of polymorphism, average numbers of alleles per locus, Shannon-Weaver information indices 
of allelic diversity, as well as values of haplotype and nucleotide diversity were similar in all regio¬ 
nal samples. But hares from both New Zealand and Britain had significantly lower genetic diversity 
than brown hares from continental Europe studied earlier. Thus, gene pool erosion likely occurred 
already in British hares, perhaps associated with their probable introduction in Roman times. The¬ 
oretically, the small number of alleles found in British brown hares could have been sampled by 
the few hares that were reported as having constituted the founder stock in New Zealand in the 
nineteenth century. As expected, rare alleles of British brown hares were absent in New Zealand. 
But drift had only a slight effect on the gene pool composition of hares in New Zealand. 
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Introduction 

In New Zealand, brown hares (Lepus euro- 
paeus ) have higher rates of ovarian tumors 
and missing posterior upper molars (M 3 ) 
than in Europe (Flux 1965, 1980; Parkes 
1988; Suchentrunk et al. 1992). This might 
result from low genetic variability as a con¬ 
sequence of a small number of founder in¬ 
dividuals (Flux 1965). Historical docu¬ 
ments suggest that brown hares released in 
New Zealand in the 19 th century by various 
Acclimatization Societies were mostly ta¬ 


ken from Phillip Island, Victoria, Australia 
(Lever 1985; Flux 1990). There, only six 
hares had built up a population of 200 indi¬ 
viduals by 1865, few years after introduc¬ 
tion (Flux 1990; Flux et al. 1990; see also 
Rolls 1969 and Lever 1985). In Australia, 
brown hares were probably first success¬ 
fully introduced in 1859 by W. Lyall on 
the shores of Western Port Bay, Victoria 
(Lever 1985), and afterwards on Phillip Is¬ 
land (Mahood 1983). In February 1864 an- 
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other nine hares were released into an en¬ 
closure near Geelong, Victoria by T. Aus¬ 
tin, who imported them from England (cf. 
Lever 1985). All Australian hares are con¬ 
sidered originating from Britain (cf., Flux 
1990), but no details as to specific regions 
are given in the available literature (Lever 
1985; Flux pers. com.). The exact numbers 
of hares that have successfully bred after 
their naturalization in Australia and New 
Zealand remain unknown. However, the 
list of importations to New Zealand pre¬ 
sented by Lever (1985) suggests limited ge¬ 
netic variability in the founder gene pool 
(see also Flux 1990). 

In this study we compared levels of genetic 
variability of hares from New Zealand and 
Britain, to test this genetic bottleneck hy¬ 
pothesis (Flux 1965). Theory and empirical 
findings (e. g., Fuerst and Maruyama 1986; 
Leberg 1992; Hartl and Pucek 1994; Tie- 
demann et al. 1997) predict a smaller effect 
of bottlenecks on multi-locus allozyme het¬ 
erozygosity than on other indicators of allo- 
zymic variability, such as the rate of poly¬ 
morphism (P), and mean number of alleles 
per locus (A). Allozyme heterozygosity 
may even increase after bottlenecks (e.g., 
Leberg 1992). Thus, we expected lower P- 
and A-values for hares from New Zealand 
than for British brown hares, whereas het¬ 
erozygosities might be similar. Particularly 
alleles with low frequencies in British 
brown hares might have not been sampled 
by the few founders in New Zealand. Drift 
effects could have caused shifts in allele fre¬ 
quencies and consequently increased genet¬ 
ic divergence between hares from Britain 
and New Zealand. We also expected a pro¬ 
nounced reduction of variability in the mi¬ 
tochondrial DNA (mtDNA), because of 
the lack of recombination in this maternally 
inhereted genome in post-bottleneck popu¬ 
lations (e. g., Giles et al. 1980; Lansman et 
al. 1981; Avise 1994; Avise and Hamrick 
1996; see e.g., Gyllensten et al. 1991 for 
paternal inheritance of mtDNA). Brown 
hares from the British Isles are convention¬ 
ally considered a separate subspecies (L. e. 
occidentalis De Winton, 1898; cf. Corbet 
and Southern 1977; Arnold 1993). They 


might be genetically somewhat distinct 
from mainland European brown hares. 
Therefore, we compared the present data 
with adjusted data sets of continental Eur¬ 
opean brown hares published earlier 
(Hartl et al.1993; Suchentrunk et al. 
2000 ). 


Material and methods 

Specimens studied 

We studied 119 hares from two regions in New 
Zealand and three in Britain (Fig. 1). In New 
Zealand, hares were collected in the Wairarapa 
region (n = 32) of the North Island in September/ 
October 1993 by J. E. C. Flux (“Landcare Re¬ 
search”, Lower Hutt), and in the Harper/Avoca 
catchment (n = 28) on the South Island in Octo¬ 
ber 1993 and March 1995 by J. Parkes (“Landcare 
Research”, Christchurch) and F. Suchentrunk. In 
Britain, collections were organized by S. Tapper 
(The Game Concervancy Trust, Fordingbridge, 
England) in February 1995 in three regions (Wilt¬ 
shire, southern England, n = 20; Loddington, Lei¬ 
cestershire, central England, n = 19; Duns, Aber¬ 
deenshire, Scotland, n = 20). Most hares were 
dissected by one of the authors (FS). They were 
sexed by inspection of the internal reproductive 
organs and aged by checking the lateral epiphy¬ 
sial protrusion of the ulna (Stroh’s sign), which 
separates young of the year (< approx. 7- 
10 months) from older hares (Suchentrunk et 
al. 1991). Liver, kidney, and spleen tissue samples 
were frozen at -20 °C until further use. 


Allozymic diversity 

We screened allelic variation at 52 hypothetical 
structural gene loci in 119 hares by standard hori¬ 
zontal starch gel electrophoresis of the following 
isozymes/-systems (isozyme/-system, abbrevia¬ 
tion, E.C. number, and corresponding structural 
gene loci in parentheses): a-glycerophosphate de¬ 
hydrogenase (GDC, 1.1.1.8, Gdc), sorbitol dehy¬ 
drogenase (SDH, 1.1.1.14, Sdh), lactate dehydro¬ 
genase (LDH, 1.1.1.27, Ldh-1,-2), malate 

dehydrogenase (MOR, 1.1.1.37, Mor-1,-2), malic 
enzyme (MOD, 1.1.1.40, Mod-1,-2), isocitrate de¬ 
hydrogenase (IDH, 1.1.1.42, Idh-1,-2), 6-phospho- 
gluconate dehydrogenase (PGD, 1.1.1.44, Pgd), 
glucose dehydrogenase (GDH, 1.1.1.47, Gdh-2), 
glucose-6-phosphate dehydrogenase (GPD, 
1.1.1.49, Gpd), glyceraldehyde-3-phosphate dehy- 
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Fig. 1 . Sampling regions of brown hares in Britain and New Zealand. Details are given in the text. 


drogenase (GAPDH, 1.2.1.12, Gapdh), xanthine 
dehydrogenase (XDH, 1.2.3.2, Xdh), glutamate 
dehydrogenase (GLUD, 1.4.1.3, Glud), NADH- 
diaphorase (DIA, 1.6.2.2., Dia-1), catalase (CAT, 
1.11.1.6, Cat), superoxide dismutase (SOD, 

1.15.1.1, Sod-1,-2), purine nucleoside phosphory- 
lase (NP, 2.4.2.1, Np), aspartate aminotransferase 
(AAT, 2.6.1.1, Aat-1,-2), hexokinase (HK, 2.7.1.1, 
Hk-1,-2,-3), pyruvate kinase (PK, 2.7.1.40, Pk-1), 
creatine kinase (CK, 2.7.3.2, Ck-1,-2), adenylate 
kinase (AK, 2.7.4.3, Ak-1,-2), phosphoglucomu- 
tase (PGM, 2.7.5.1, Pgm-1,-2,-3), esterases (ES, 

3.1.1.1, Es-1; 4.2.1.1, Es-D), acid phosphatase 
(ACP, 3.1.3.2, Acp-1,-2,-3), fructose-1,6-dipho¬ 
sphatase (FDP, 3.1.3.11, Fdp-1), peptidases (PEP, 
3.4.11, Pep-1,-2), guanine deaminase (GDA, 

3.5.4.3, Gda), adenosine deaminase (ADA, 

3.5.4.4, Ada-2,-3), aldolase (ALDO, 4.1.2.13, 
Aldo), fumarate hydratase (FH, 4.2.1.2, Fh), aco- 
nitase (ACO, 4.2.1.3, Aco-1,-2), mannose phos¬ 
phate isomerase (MPI, 5.3.1.8, Mpi), glucose 
phosphate isomerase (GPI, 5.3.1.9, Gpi-1,-2). This 
suite of isozyme loci is identical to the loci 
screened by Hartl et al. (1993) for 20 regional 
samples of brown hares from central Europe, ex¬ 
cept for Acy-1 and /?-Gal, which were not 
screened presently. 


In tissue preparation, electrophoresis, and protein 
specific staining we followed Grillitsch et al. 
(1992). For resolving allelic variants we made di¬ 
rect side-by-side comparisons of migrating allo- 
zymes on the same gels by including samples of 
brown hares studied earlier in our laboratory, 
and adopted the allele nomenclature of Hartl et 
al. (1993) and Suchentrunk (1993). We deter¬ 
mined genotypes at polymorphic loci from zymo¬ 
grams according to principles of allozyme electro¬ 
phoresis (Harris and Hopkinson 1976; Rothe 
1994). In some individuals, however, we could 
not genotype all polymorphic loci due to poor re¬ 
solution. 

We used the BIOSYS-1 pc package, release 1.7 
(Swofford and Selander 1989) to calculate al¬ 
lele frequencies, average heterozygosity (H 0 - ob¬ 
served, H e - expected), proportion of poly¬ 
morphic loci (P, 99% criterion), mean number of 
alleles per locus based on all 52 loci (A), and for 
exact tests of deviations of observed genotypes at 
polymorphic loci from Hardy-Weinberg expecta¬ 
tions. As additional index of genetic diversity, we 
calculated the Shannon-Weaver information in¬ 
dex (H'; see Hedrick 1985) for each regional sam¬ 
ple (RS) as the sum of locus-specific information 
indices based on allele frequencies. 
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We tested variation of locus-specific allele fre¬ 
quencies between pairs of regional samples by ex¬ 
act Fisher’s tests and based significance decisions 
on Sequential Bonferroni procedure (Rice 1989), 
to account for multiple testing. Also, we based ex¬ 
act Fisher’s tests of variation of allele frequencies 
between the two age classes (young of the year 
vs. older animals) and sex on Sequential Bonfer¬ 
roni procedure. In order to evaluate possible drift 
effects by e. g., founder events and/or long-term 
low effective population sizes, we calculated pair¬ 
wise genetic D values (Nei 1978), modified Ro¬ 
gers’ distances, and fixation indices (F ST ) (Wright 
1978) among RSs by the BIOSYS-1 pc program 
package. Furthermore, we calculated RS-specific 
Fis -values with this software, to check whether 
or not possible low genetic variability might result 
from regional inbreeding. For direct comparison 
with allozymic variability of brown hares from 
20 local samples from Austria (Hartl et al.1993) 
and eight regional samples from Bulgaria (Su- 
chentrunk et al. 2000), we adjusted all samples 
to 49 loci by omitting the Acy-1, /?-Gal, Ada-2, 
Ada-3, and Dia-1 loci. All these data sets were 
produced in our laboratory by using respective 
marker individuals on the gels, which enabled a 
direct comparison. 


MtDNA-RFLP analysis 

We used liver tissue samples of 36 hares (NZ- 
W = 10, NZ-H = 7, UK-W = 3, UK-L = 13, UK- 
D = 3) to isolate mtDNA by CsCl/EtBr density 
gradient ultra-centrifugation and removing EtBr 
by isoamyl alcohol extraction, basically following 
the protocols of Lansman et al. (1981) and Rick- 
wood (1987); for details see Hartl et al. (1993) 
and Nadlinger (1994). We digested total mtDNA 
of each hare with the following set of 16 hexanu- 
cleotid-recognizing type II restriction endonu¬ 
cleases (Roche): Apal, BamHI, Bell, Bglll, Clal, 
Dral, EcoRI, EcoRV, Hindlll, Hpal, PstI, PvuII, 
SacI, Xbal, Xhol, and Xmnl. The resultant frag¬ 
ments were separated electrophoretically (80 V, 
two hours) in 0.7% agarose gels with 0.5 pg 
EtBr/ml and visualized in UV light. In order to 
validate length estimations of > 4 kb long frag¬ 
ments, electrophoresis was continued for two ad¬ 
ditional hours and fragment measurements were 
repeated. Fragment length variants < 0.4 kb could 
not be detected by our procedure. Fragment 
lengths were determined by comparing with 
Lambda phage DNA digested with HindiII and 
a 100 bp-ladder (1500-500 bp range). We com¬ 
pared fragments and cleaveage sites deduced 
from enzyme-specific fragment patterns with 


those already found in brown hares from Austria 
(Hartl et al. 1993; Nadlinger 1994). Because 
our enzyme set was identical with that used al¬ 
ready for Austrian brown hares, we could com¬ 
pare our restriction morphs, haplotypes, and in¬ 
dices of mtDNA variability directly with brown 
hare mtDNA data published by these authors 
(Hartl et al. 1993). We calculated haplotype di¬ 
versity (h) and nucleotide diversity ( n ) (Nei and 
Li 1979; Nei 1987) to describe RS-specific gene 
pool variability. We calculated pairwise net nu¬ 
cleotide diversity between RSs based on cleave¬ 
age site variations (e. g., Avise 1994) to obtain es¬ 
timates of mtDNA differentiation among RSs. 
We compared frequencies of the standard haplo¬ 
type (i. e., the by far most common type I; Hartl 
et al. 1993) and other haplotypes (i. e., all others 
aggregated) between NZ and UK hares by a 
one-tailed exact Fisher test (hypothesizing lower 
variability in NZ than UK hares). Finally, we 
compared RS-specific h- and ^-values of NZ and 
UK hares with the respective values of 18 Aus¬ 
trian brown hare samples (Hartl et al. 1993) by 
Mann-Whitney U-tests, basing significance deci¬ 
sions on Sequential Bonferroni procedure. 


Results 

We found di-allelic variation at six loci 
(Tab. 1) and a significant excess of homozy¬ 
gotes at the Es-1 locus in the NZ-W sample. 
Except for Acp-1 50 and Mpi 84 , all alleles 
were found earlier in brown hares from 
central Europe. Allele frequencies did not 
depend on age category or sex. The indices 
of genetic variability (H 0 , H e , P, A, H') are 
listed in table 1. Overall and locus-specific 
inbreeding coefficients (F IS ) are given in ta¬ 
ble 2, separately for each RS. Nei’s (1978) 
unbiased D values, Rogers’ modified dis¬ 
tances, and the fixation coefficients (Fst) 
appear in table 3, along with the signifi¬ 
cances of pairwise comparisons of allele fre¬ 
quencies at one or more loci. In table 4, the 
indices of allozymic variability (based on 
49 loci) of the UK and NZ samples are 
compared to those of 20 Austrian and eight 
Bulgarian regional samples of brown hares 
studied earlier (Hartl et al.1993, Suchen- 
trunk et al. 2000). 

Thirtyfive (97.2%) hares had the standard 
mtDNA haplotype I, that was already 
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Table 1. Allele frequencies at polymorphic loci, sample-specific indices of genetic variability, and mean inbreed¬ 
ing coefficients (F IS ) in brown hares from New Zealand (NZ) and Britain (UK); for acronyms see figure 1. Allele 
designations correspond to those given in Hartl et al. (1989, 1993), Suchentrunk (1993), Suchentrunk et al. 
(1998), and Suchentrunk et al. (2000). n = number of individuals screened for each locus and regional sample. In¬ 
dices of gene pool variability (based on 52 loci): H 0 = observed average heterozygosity, H e = expected average 
heterozygosity, P( 99 ° /o ) = rate of polymorphism (99% criterion), A = average number of alleles per locus, 
H' = Shannon Weaver information index. 


Locus 

Allele 

NZ-W 


NZ-H 

UK-W 

UK-L 

UK-D 

Idh-2 

n 

32 


28 

20 

19 

20 


100 

1.0 


1.0 

0.975 

1.0 

1.0 


130 

0.0 


0.0 

0.025 

0.0 

0.0 

Pep-2 

n 

32 


28 

20 

19 

20 


100 

0.781 


0.536 

0.850 

1.0 

0.875 


104 

0.219 


0.464 

0.150 

0.0 

0.125 

Acp-1 

n 

32 


28 

20 

19 

20 


100 

1.0 


1.0 

1.0 

1.0 

0.975 


50 

0.0 


0.0 

0.0 

0.0 

0.025 

Ada-2 

n 

32 


6 

18 

17 

18 


100 

0.813 


0.750 

0.694 

0.588 

0.639 


75 

0.187 


0.250 

0.306 

0.412 

0.361 

Es-1 

n 

32 


28 

20 

19 

20 


-100 

0.781 


1.0 

0.600 

0.632 

0.650 


-75 

0.219 


0.0 

0.400 

0.368 

0.350 

Mpi 

n 

32 


28 

20 

19 

20 


100 

0.938 


0.964 

1.0 

1.0 

1.0 


84 

0.062 


0.036 

0.0 

0.0 

0.0 

H 0 


0.019 


0.013 

0.026 

0.017 

0.023 

He 


0.022 


0.019 

0.024 

0.019 

0.023 

P ( 99 %) 


7.69 


5.77 

7.69 

3.85 

7.69 

A 


1.08 


1.06 

1.08 

1.04 

1.08 

H' 


1.767 


1.408 

1.829 

1.336 

1.795 

Table 2. Locus-specifc (unbiased) heterozygosities in 

% (upper values) and inbreeding coefficients (lower va- 

lues) as well as overall inbreeding coefficients (F IS ) for regional samples of brown hares from New Zealand (NZ) 

and Britain (UK). 








NZ-W 


NZ-H 


UK-W 

UK-L 

UK-D 

Idh-2 

_ 


- 


5.0 

_ 

- 






-0.026 



Pep-2 

34.7 


50.6 


26.2 

- 

22.4 


+0.086 


+0.138 

-0.176 


-0.143 

Acp-1 

- 


- 


- 

- 

5.0 








-0.026 

Ada-2 

31.0 


40.7 


43.7 

49.9 

47.5 


-0.231 


+0.556 

-0.178 

+0.029 

-0.084 

Es-1 

34.7 


- 


49.2 

47.8 

46.7 


+0.0451 



-0.042 

+0.095 

+0.121 

Mpi 

11.9 


7.0 


- 

- 

- 


-0.067 


-0.067 




overall 








Fis 

+0.095 


+0.292 

-0.117 

+0.061 

-0.014 
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found in the majority of central European 
brown hares (EIartl et al. 1993), and only 
one hare of the NZ-H sample had a new 
haplotype. This haplotype deviated from 
haplotype I by only one additional cleave- 
age site, produced by Xbal at position 14.8 
kb of the restriction map published in 
Hartl, et al. (1993). Values of haplotype 
(h) and nucleotide diversity ( n ) were zero 
for all RSs except for the NZ-H sample; 


that had an h-value of 28.57% and a 7r-value 
of 0.049%. Haplotype frequencies did not 
differ significantly between the RSs from 
Britain and New Zealand. Also, RS-specific 
h- and ^-values of the samples from Britain 
and New Zealand were not significantly 
lower (p>0.05, one-tailed Mann-Whitney 
U-tests) than in brown hares from 20 Aus¬ 
trian localities (cf. Hartl et al. 1993). Va¬ 
lues of pairwise net nucleotide diversities 


Table 3. Matrix of pairwise genetic distances and fixation indices (F ST ) among regional samples of brown hares 
from New Zealand and Britain. Nei's (1978) unbiased D values (first row) and modified Rogers' distances (Wright 
1978) (second row) above the diagonal and F S t values below. F ST values were considered differing significantly 
from zero with significant allele frequencies for at least one locus in pairwise comparisons (significance based on 
Sequential Bonferroni procedure; nominal a = 0.05, 45 tests). Significance (sig) or no significance (n. s.) is indi¬ 
cated below each F ST value. For acronyms of regional samples see figure 1. 


regional 

samples 

(i) 

(2) 

(3) 

(4) 

(5) 

NZ-W (1) 

_ 

0.001 

0.001 

0.001 

0.001 



0.047 

0.033 

0.049 

0.034 

NZ-H (2) 

0.052 

- 

0.004 

0.007 

0.004 


sig 


0.071 

0.085 

0.070 

UK-W (3) 

0.024 

0.109 

- 

0.000 

0.000 


n. s. 

sig 


0.026 

0.019 

UK-L (4) 

0.057 

0.167 

0.016 

- 

0.000 


sig 

sig 

n. s. 


0.019 

UK-D (5) 

0.026 

0.106 

0.003 

0.009 

- 


n. s. 

sig 

n. s. 

n. s. 



Table 4. Comparison of indices of allozymic and mtDNA variability of regional samples of brown hares from New 
Zealand, Britain, Austria, and Bulgaria. Allozyme data are based on 49 loci and mtDNA RFLP data on 16 restric¬ 
tion endonucleases (see Material and methods). Means and range (in parentheses) of observed (H 0 ) and expected 
(H e ) average heterozygosity, Shannon-Weaver diversity index (H'), rate of polymorphism (P - 99% criterion), 
average number of alleles per locus (A), haplotype (h) and nucleotide diversity (n) are given for each group of 
local samples, sig = significance as determined by Mann Whitney tests (d. f. = 1) and sequential Bonferroni pro¬ 
cedure (a = 0.05), n. s.= not significant. 


index 

2 NZ and 3 UK regional 
samples (this study) 

20 Austrian (Hartl et al. 1993) and 

8 Bulgarian (Suchentrunk et al. 

2000) regional samples 

sig. 

H 0 (%) 

1.26 (0.9-1.7) 

2.83 (1.7-4.2) 

<0.0001 

He (%) 

1.4 (1.0-1.7) 

3.06 (2.3-4.7) 

<0.0001 

H' 

1.029 (0.66-1.29) 

2.431 (1.73-3.66) 

<0.0001 

P (99%) 

4.90 (2.04-6.12) 

11.15 (8.16-16.33) 

<0.0001 

A 

1.048 (1.02-1.06) 

1.144 (1.08-1.2) 

<0.0001 

He/P 

0.314 (0.245-0.490) 

0.279 (0.172-0.451) 

n.s. 

h (%) 

5.7 (0-28.57) 

15.9 (0-69.9) 1 

n.s. 

n (%) 

0.0098 (0-0.049) 

0.0351 (0-0.184) 1 

n. s. 


1 mtDNA values calculated only from data of 18 Austrian regional samples 
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were zero for the British and NZ-W sam¬ 
ples, and amounted to 0.0041% for all pairs 
involving the NZ-H sample. Respective 
pairwise values for the presently studied 
UK and NZ hares and the earlier studied 
brown hares from 20 Austrian localities 
(Hartl et al. 1993) ranged between 0.0- 
0.093%; this was within the range (0.0- 
0.113%) of net nucleotide diversity be¬ 
tween local samples of brown hares from 
Austria (calculated from data produced in 
our laboratory, see Hartl et al. 1993). 


Discussion 

The levels of cross nuclear and mtDNA 
variability of brown hares from New Zea¬ 
land and Britain are similar. This is in con¬ 
tradiction to the hypothesis of reduced gene 
pool variability in hares from New Zealand, 
owing to an assumed bottleneck during the 
period of their introduction. Possibly, un¬ 
recorded liberations of hares from diverse 
provenences in Europe in addition to those 
from the reportedly small number of foun¬ 
ders in Australia (Lever 1985) have in¬ 
creased the effective population size during 
the founder period in New Zealand. But 
such additional imports from Europe to 
New Zealand do not seem very likely, given 
the acclimatized hares already available in 
Australia for the Acclimatization Compa¬ 
nies. 

The essential point to explain the absence 
of reduced genetic variability in hares from 
New Zealand as compared to British hares 
is that brown hares from both New Zealand 
and Britain exhibit clearly lower allozymic 
variability than continental European po¬ 
pulations. Theoretically, the small number 
of allozymic alleles found presently in the 
British brown hares could have been 
sampled by only few individuals. All com¬ 
mon alleles (i. e., those with relative electro¬ 
phoretic mobility 100/-100) of the British 
samples were also common in brown hares 
from diverse regions of continental Europe 
(Hartl et al. 1989, 1990, 1992, 1993, 1994; 
Suchentrunk et al. 2000; Suchentrunk et 
al. unpubl. data). But only 16% of all var¬ 


iant alleles found so far in continental Eur¬ 
opean brown hares with the same set of loci 
(cf. Hartl et al. 1994; Suchentrunk et al. 
2000) occurred in the British samples. Quite 
several with wide distribution in continental 
Europe (Pgd 129 , Hk-2 67 , Es-1^ 42 , Es-L 108 , 
Es-D 141 , Mp 126 ) are likely absent in British 
brown hares. And only one (7.1%) of all 
mtDNA haplotypes found so far in brown 
hares from central and southeastern Eur¬ 
ope (Hartl et al. 1993, 1994; Nadlinger 
1994; Suchentrunk, unpubl. data) could 
be detected in British brown hares. But 
quite a number of regional samples of Aus¬ 
trian brown hares did also not display any 
mtDNA variability (Hartl et al. 1993). 
British sample sizes are probably too low 
for such a comparison. Nevertheless, ab¬ 
sence of any mtDNA variability of the pre¬ 
sently studied brown hares from Britain 
and the fact that only the standard Eur¬ 
opean mtDNA haplotype I was found in 
Britain, as opposed to its significantly 
(p< 0.0001, exact Fisher’s test) lower fre¬ 
quency on the continent, agrees with the in¬ 
terpretation of generally reduced genetic 
variability in British brown hares. 

The low genetic variability of British brown 
hares might result from an “ancient" popu¬ 
lation bottleneck or long-term low effective 
population sizes associated with their colo¬ 
nization history. It has been hypothesized 
that brown hares were introduced in Ro¬ 
man times (cf. Corbet 1986; Arnold 
1993). Deliberate releases of only few indi¬ 
viduals or occasional escapes from farms in 
Roman times or earlier during the Meso¬ 
lithic (?) or Neolithic occupation periods 
(cf. Jones and Keen 1993) could have re¬ 
sulted in a poor gene pool variability of the 
pioneer population. In addition, a low sur¬ 
vival rate owing to a high predation pres¬ 
sure by foxes and other predators, patho- 
genes, adverse weather conditions etc., 
could have hampered a quick population 
growth in the wild. Long-term low effective 
population size effectively reduces allelic 
variabililty and heterozygosity (e. g., Chak- 
raborty and Nei 1977; Hedrick 1985). Con¬ 
trary to the hypothesis of a deliberate intro¬ 
duction of brown hares to Britain in Roman 
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times, Yalden (1982) and Roberts (1994) 
list brown hares along with mountain hares 
(Leptus timidus) as native to Britain since 
late-glacial times (see also Grigson 1983 
for the Later Mesolithic excavation site of 
Cherhill, Jones and Keen 1993, and cita¬ 
tions therein). In this case, brown hares 
could have lost genetic diversity at low den¬ 
sities under adverse late-glacial climate, or 
later on in small isolated pockets. Stuart 
(1982), however, lists “Lepus sp.” as an ele¬ 
ment of the early Flandrian fauna of Star 
Charr, and considers only L. timidus defi- 
nitly recorded as fossil in Britain (cf., May- 
hew 1975). As to our knowledge, fossil evi¬ 
dence of Lepus europaeus is uncertain in 
north-central and north-western mainland- 
Europe during the late-glacial period and 
early Flandrian, before the formation of 
the Channel (Strait of Dover, c.8000 years 
BP cf. e. g., Jones and Keen 1993). Brown 
hares might not have roamed mid-latitude 
Europe during the early Flandrian and not 
have managed to arrive there before the 
formation of the Channel (see also Corbet 
1986). Alternatively, genetic diversity could 
have been reduced in fragmented popula¬ 
tions with long-term low densities in subop- 
timal habitats after woodland regeneration 
in the post-Roman period. But it seems that 
the distribution of wooded and non-wooded 
land in Britain has not been altered much in 
the post-Roman period (Roberts 1994). 

The low allozymic diversity of British 
brown hares unlikely results from recent in- 
breeding in regional populations. This is in¬ 
dicated by quite normal locus-specific he¬ 
terozygosities, lack of heterozygote defi¬ 
ciencies, and the generally low or even ne¬ 
gative inbreeding coefficients. This inter¬ 
pretation is also supported by the absence 
of a significant gene pool substructuring. 
Allele frequencies of British brown hares 
do not differ much across regions. Insignifi¬ 
cant fixation indices (F S t) and genetic dis¬ 
tance values indicate absence of drift effects 
among local populations. Obviously, the 
low amount of genetic variability contained 
in British brown hares is partitioned among 
individuals within local populations rather 
than across larger geographic ranges. The 


Acp-1 50 allele occurred exclusively in the 
Scottish sample and was not found in a 
large number of brown hares from many re¬ 
gions in continental Europe. It might result 
from a recent local mutation. 
Differentiation between the gene pools of 
British and central European (Austrian) 
brown hares (cf., Hartl et al. 1993) is negli¬ 
gible, and virtually nil when based on the 
detected mtDNA haplotypes. The corre¬ 
sponding pairwise Nei’s (1978) unbiased D 
values (0.000-0.005) fall within the range 
encountered among local populations from 
Austria (e. g., Hartl et al. 1993). Despite 
their conventional subspecific position 
(L. e. occidentals, De Winton, 1898), Brit¬ 
ish brown hares represent only a genetically 
depauperate version of brown hares from 
continental Europe. This is confirmed by 
the absence of any other mtDNA haplotype 
apart from haplotype I, which has a phylo- 
genetically central position in brown hares 
and is most widespread in central Europe 
(Hartl et al. 1993; Suchentrunk, unpubl. 
data). As indicated by the low genetic dis¬ 
tances, no significant genetic drift has oc¬ 
curred between British and continental 
European brown hare populations. As re¬ 
gards mtDNA, we cannot draw any conclu¬ 
sions on differentiation among the regions 
studied in Britain, because of too low sam¬ 
ple sizes for two regions. Also, the present 
data do not permit any conclusion as to ori¬ 
gins of the British brown hares. This might 
be achieved by comparing regional popula¬ 
tion samples and continental European 
samples with a highly resolving molecular 
marker system (e. g., microsatellites, mito¬ 
chondrial d-loop sequences). 

The lack of distinct gene pool structuring 
and the poor allozymic diversity of British 
brown hares suggest that a relatively high 
proportion of variant alleles was sampled 
by the few individuals that supposedly con¬ 
stituted the founder populations in New 
Zealand. Theoretically, all currently found 
alleles of the British brown hares could 
have been contained in the founders. 

In New Zealand, hares spread rapidly after 
their introduction (Flux 1990). Presently, 
they occupy large ranges including practi- 
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cally all types of habitats on the North and 
South Islands. Although they commonly do 
not reach those high densities as in Europe, 
they even have established successful popu¬ 
lations in sub-alpine (Parkes 1984) and 
rather harsh alpine environments (Flux 
1967). Their successful spread in New Zea¬ 
land might have been particularly favoured 
by a low level of parasitic or infectuous dis¬ 
eases, a long-term low level of agrochemis¬ 
try, and relatively little intensification of 
farming machinery, little road mortality, 
low predation pressure on adults etc. (see 
also Flux 1990). This rapid increase of the 
founder populations on both main islands, 
together with repeated imports and subse¬ 
quent translocations within and between 
the islands by diverse naturalization compa¬ 
nies (Flux 1990) probably prevented severe 
loss of genetic variability. Generally, quick 
population increases following a bottleneck, 
as was reported for the founder populations 
of Phillip Island (Rolls 1969) and New 
Zealand (Flux 1990), are counteracting re¬ 
duction of allelic variability in founder po¬ 
pulations (e.g., Nei et al. 1975). 

As expected, particularly rare alleles (Idh- 
2 130 , Acp-1 50 ) of British brown hares are 
absent in hares from New Zealand, whereas 
the common British alleles are common in 
New Zealand too. Interestingly, the Mpi 84 
allele of both populations from New Zea¬ 
land was neither present in the British 
brown hares nor in any of over 900 brown 
hares from diverse regions of continental 
Europe studied so far (Hartl et al. 1994; 
Suchentrunk et al. 2000; Suchentrunk, 
unpubl. data). It might stem from a recent 
mutation in New Zealand or Australia. Al¬ 
ternatively, it could be present in regions 
of Britain or anywhere else in Europe that 
were not yet sampled, but wherefrom hares 
were shiped to Australia. The same inter¬ 
pretations apply to the exclusive occurrence 
of mtDNA haplotype V in the Harpa/Avoca 
catchment on the South Island of New 
Zealand. This haplotype is phylogenetically 
closely related to the basal mtDNA haplo¬ 
type I of European brown hares, and its 
evolution can be explained by only one 
base substitution. 


The two local populations on both main is¬ 
lands of New Zealand have differentiated 
little since their foundations over one hun¬ 
dred years ago. The nuclear and mitochon¬ 
drial gene pools of hares from New Zealand 
are very similar to those of British brown 
hares. Lack of data of allele frequencies in 
the source population, unknown effective 
population sizes, and repeated transloca¬ 
tions of hares within and between the North 
and South Islands, that might have changed 
allele frequencies, prevented us to compare 
the presently observed allele frequencies 
with theoretical frequencies as resulting 
from drift simulations (see e.g., FitzSim- 
mons et al. 1997). But both Nei’s (1978) D 
and FsT-values between British and New 
Zealand populations largely fall within the 
ranges encountered among regional popula¬ 
tions in Europe (D = 0.000-0.019, F S t = 
0.00-0.124, re-calculated from adjusted data 
sets of 49 loci, cf. Hartl et al. 1993; Suchen¬ 
trunk et al. 2000). In fact, only the F S i-va- 
lue for the NZ-H and UK-L samples slightly 
exceeds the range for regional samples from 
mainland Europe. This marginally increased 
level of differentiation between these two 
populations apparently results from the ab¬ 
sence of the Pep-2 104 allele in UK-L and 
the Es-1 -75 allele in NZ-H. But in general, 
we cannot see any marked reorganization 
in the cross gene pool of hares from New 
Zealand as compared to British brown 
hares. Rapid population increases in the lo¬ 
cal founder populations, repeated releases, 
and multiple transfers of hares by naturali¬ 
zation companies (Flux 1990) likely pre¬ 
vented strong genetic drift. However, in 
spite of apparently well adapted populations 
in Britain and New Zealand (Suchentrunk 
et al. 1998), brown hares from these coun¬ 
tries harbour less genetic resources for 
adaptation to future environments than po¬ 
pulations from mainland Europe. 
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Zusammenfassung 

Geringe ALlozym- und mtDNA-Variabilitat bei neuseelandischen und britischen Feldhasen 
(Lepus europaeus) - eine Folge von Flaschenhalsen? 

Untersucht wurde die Allozym- und mtDNA-Variabilitat bei britischen und neuseelandischen Feldha¬ 
sen (Lepus europaeus ), um die Hypothese eines Variabilitatsverlustes bei neuseelandischen Hasen 
infolge eines Flaschenhalsereignisses („Griinder-Effekt") bei ihrer Einburgerung im 19. Jahrhundert 
zu iiberprufen. Bei 119 Hasen aus drei britischen und zwei neuseelandischen Stichprobengebieten 
wurde mittels horizontaler Starkegelelektrophorese ihre allelische Variability an 52 Strukturgenloci 
ermittelt. Bei 36 Hasen wurde anhand von 16 6-Basenschneidenden Endonukleasen der Restrik- 
tionsfragmentlangenpolymorphismus (RFLP) in der gesamten mitochondrialen DNA (mtDNA) analy- 
siert. Alle Indices der Allozym- und mtDNA-Variabilitat lagen bei den neuseelandischen Stichproben 
im Bereich jener der britischen Feldhasen. Somit kann die Hypothese zu Verlusten an genetischer 
Variability bei neuseelandischen Hasen im Vergleich zu den britischen nicht aufrecht erhalten wer- 
den. Ausgepragte genetische Drifteffekte zwischen britischen und neuseelandischen Feldhasen 
konnten ebenfalls nicht festgestellt werden. Jedoch zeigten alle fiinf untersuchten regionalen 
Stichproben significant geringere Allozym-Variabilitat als Feldhasen vom europaischen Kontinent. 
Dies laftt die Vermutung eines Verlustes der genetischen Variability bei den britischen Feldhasen 
im Verlaufe ihrer Besiedlungsgeschichte zu. Die nominell als eigene Unterart gefuhrten britischen 
Feldhasen (L. e. occidentals, De Winton; 1898) erwiesen sich lediglich als genetisch verarmte Var- 
ianten der kontinental-europaischen Feldhasen. 

References 

Arnold, H. R. (1993): Atlas of mammals in Brit¬ 
ain. ITE research publ. 6. London: HMSO. 

Avise, J. C. (1994): Molecular Markers. Natural 
History and Evolution. New York: Chapman 
and Hall. 

Avise, J. C.; Hamrick, J. C. (1996): Conservation 
Genetics. Case Histories from Nature. New 
York: Chapman and Hall. 

Chakraborty, R.; Nei, M. (1977): Bottleneck ef¬ 
fects on average heterozygosity and genetic 
distance with the stepwise mutation model. 

Evolution 31, 347-356. 

Corbet, G. B. (1986): Relationships and origins of 
the European lagomorphs. Mammal Rev. 16, 

105-110. 


Corbet, G. B.; Southern, H. N. (1977): The 
Handbook of British Mammals. Sec. Ed.; Ox¬ 
ford: Blackwell. 

FitzSimmons, N.; Buskirk, S. W.; Smith, M. H. 
(1997): Genetic changes in reintroduced 
Rocky Mountain Bighorn Sheep populations. 
J. Wildl. Managem. 61, 863-872. 

Flux, J. E. C. (1965): Incidence of ovarian tumors 
in hares in New Zealand. J. Wildl. Managem. 
29, 622-624. 

Flux, J. E. C. (1967): Hare numbers and diet in 
an alpine basin in New Zealand. Proc. New 
Zeal. Ecol. Soc. 14, 27-33. 

Flux, J. E. C. (1980): High incidence of missing 
posterior upper molars in hares (Lepus euro- 


58 


F. SUCHENTRUNKetal. 


paeus) in New Zealand. New Zeal. J. Zool. 7, 
257-259. 

Flux, J. E. C. (1990): Brown hare. In: The Hand¬ 
book of New Zealand Mammals. Ed. by C. M. 
King. Auckland: Oxford University Press. Pp. 
162-172. 

Flux, J. E. C.; Duthie, A. G.; Robinson, T. J.; Chap¬ 
man, J. A. (1990): Exotic populations. In: Rab¬ 
bits, Hares and Pikas. Status Survey and Conser¬ 
vation Action Plan. ED. by J. A. Chapman and 
J. E. C. Flux. Gland: IUCN. Pp. 147-153. 

Fuerst, P. A.; Maruyama, T. (1986): Considera¬ 
tions on the conservation of alleles and of 
genic heterozygosity in small managed popu¬ 
lations. Zoo Biology 5,171-179. 

Giles, R E.; Blanc, H.; Cann, H. M.; Walla¬ 
ce, D. C. (1980): Maternal inheritance of hu¬ 
man mitochondrial DNA. Proc. Natl. Acad. 
Sci., USA 77, 6715-6719. 

Grigson, C. (1983): Mesolithic and neolithic ani¬ 
mal bones. In: Excavations at Cherhill, North 
Wiltshire, 1967. Proceed. Prehistoric Soc. 49. 
Ed. by J. G. Evans, and I. F. Smith. Pp. 64-72. 

Grillitsch, M.; Hartl, G. B.; Suchentrunk, F.; 
Willing, R. (1992): Allozyme evolution and 
the molecular clock in the Lagomorpha. Acta 
Theriol. 37,1-13. 

Gyllensten, U.; Wharton, U. D.; Josefson, A.; 
Wilson, A. C. (1991): Paternal inheritance of 
mitochondrial DNA. Nature 352, 255-257. 

Harris, H.; Hopkinson, D. A. (1976): Handbook 
of Enzyme Electrophoresis in Human Genet¬ 
ics. Amsterdam: North-Holland Publ. Co. 

Hartl, G. B.; Pucek, Z. (1994): Genetic depletion 
in the European bison {Bison bonasus) and 
the significance of electrophoretic heterozyg¬ 
osity for conservation. Conserv. Biol. 8, 167- 
174. 

Hartl, G. B.; Willing, R.; Nadlinger, K. (1994): 
Allozymes in mammalian population genetics 
and systematics: Indicative function of a mar¬ 
ker system reconsidered. In: Molecular Ecol¬ 
ogy and Evolution: Approaches and Applica¬ 
tions. Ed. by B. Schierwater, B. Streit, 
G. P. Wagner, and R. DeSalle. Basel: Bir- 
khauser. P. 299-310. 

Hartl, G. B.; Suchentrunk, F; Nadlinger, K.; 
Willing, R. (1993): An integrative analysis of 
genetic differentiation in the brown hare Le- 
pus europaeus based on morphology, allo¬ 
zymes and mitochondrial DNA. Acta Theriol. 
38, 33-57. 

Hartl, G. B.; Suchentrunk, F.; Willing, R.; 
Grillitsch, M. (1989): Biochemisch-genet- 
ische Variabilitat und Differenzierung beim 
Feldhasen ( Lepus europaeus) in Niederoster- 
reich. Wien. Tierarztl. Mschr. 76, 279-284. 


Hartl, G. B.; Markowski, J.; Kovacs, G.; Gril¬ 
litsch, M.; Willing, R (1990): Biochemical var¬ 
iation and differentiation in the brown hare {Le¬ 
pus europaeus) of Central Europe. Z. 
Saugetierkunde 55,186-193. 

Hartl, G. B.; Markowski, J.; Swiatecki, A.; Ja- 
niszewski, T.; Willing, R. (1992): Genetic di¬ 
versity in the Polish brown hare Lepus euro¬ 
paeus Pallas, 1778: implications for 
conservation and management. Acta Theriol. 
37, 15-25. 

Hedrick, P. W. (1985): Genetics of Populations. 
Boston: Jones and Bartlett Publ., Inc. 

Jones, R. L.; Keen, D. H. (1993): Pleistocene En¬ 
vironments in the British Isles. London: Chap¬ 
man and Hall. 

Lansman, R. A.; Shade, R. O.; Shapira, J. F.; Avi- 
se, J. C. (1981): The use of restriction endonu¬ 
cleases to measure mitochondrial DNA se¬ 
quence relatedness in natural populations. Ill 
Techniques and potential applications. J. Mol. 
Evol. 17, 214-226. 

Leberg, P. L. (1992): Effects of population bottle¬ 
necks on genetic diversity as measured by al¬ 
lozyme electrophoresis. Evolution 46, 477- 
494. 

Lever, C. (1985): Naturalized Mammals of the 
World. London, New York: Longman. 

Mahood, I. T. (1983): Brown hare {Lepus capen- 
sis). In: The Australian Museum Complete 
Book of Australian Mammals. Ed. by R. Stra- 
han. London: Angus and Robertson Publ. 
P. 480. 

Mayhew, D. F. (1975): The Quaternary history of 
some British Rodents and Lagomorphs. Ph. 
D. Thesis, University of Cambridge. 

Nadlinger, K. F. (1994): Restriktionsfragmen- 
tlangenpolymorphismus mitochondrialer 
DNA in Lepus europaeus, Cervus elaphus und 
Bison bonasus. PhD Diss., University of Vien¬ 
na. 

Nei, M. (1978): Estimation of average heterozyg¬ 
osity and genetic distance from a small num¬ 
ber of individuals. Genetics 89, 583-590. 

Nei, M. (1987): Molecular Evolutionary Genetics. 
New York: Columbia Univ. Press. 

Nei, M.; Li, W. H. (1979): Mathematical model 
for studying genetic variation in terms of re¬ 
striction endonucleases. Proc. Nat. Acad. Sci., 
U.S.A. 76, 5269-5273. 

Nei, M.; Maruyama, T.; Cakraaborty, R. (1975): 
The bottleneck effect and genetic variability 
in populations. Evolution 29,1-10. 

Parkes, J. P. (1984): Home ranges of radio-tele- 
metred hares {Lepus capensis) in a sub-alpine 
population in New Zealand: implications for 
control. Acta Zool. Fenn. 171, 279-281. 


Little allozyme and mtDNA variability in brown hares 


59 


Parkes, J. P. (1988): Missing posterior upper mo¬ 
lars in Lepus europaeus in New Zealand. La- 
gomorph Newsl. 7,12. 

Rickwood, D. (1987): Centrifugation. A Practical 
Approach. Sec. Ed. Oxford: IRL Press. 

Rice, W. S. (1989): Analyzing tables of statistical 
tests. Evolution 43, 223-225. 

Roberts, N. ( 1994): The Holocene. An Environ¬ 
mental History. Oxford: Blackwell. 

Rolls, E. C. (1969): They All Run Wild. The 
Story of Pests on the Land of Australia. Syd¬ 
ney: Angus and Robertson. 

Rothe, G. M. (1994): Electrophoresis of Enzymes. 
Laboratory Methods. Springer LAB Manual. 
Berlin: Springer. 

Stuart, A. J. (1982): Pleistocene Vertebrates in 
the British Isles. London, New York: Long¬ 
man. 

Suchentrunk, F. (1993): Variability of minor 
tooth traits and allozymic diversity in brown 
hare Lepus europaeus populations. Acta 
Theriol. 38, Suppl., 2, 59-69. 

Suchentrunk, F.; Willing, R.; Hartl, G. B. 
(1991): On eye lens weights and other age cri¬ 
teria of the Brown hare ( Lepus europaeus 
Pallas, 1778). Z. Saugetierkunde 56, 365-374. 

Suchentrunk, F.; Markowski, J.; Janiszewski, T.; 
Hartl, G. B. (1992): Dental and cranial 
anomalies in Austrian and Polish brown hares 
Lepus europaeus populations. Acta Theriol. 
37,241-257. 

Suchentrunk, F; Michailov, G; Markov, G.; 
Haiden, A. (2000): Population genetics of 
Bulgarian brown hares Lepus europaeus : allo¬ 


zymic diversity at zoogeographical crossroads. 
Acta Theriol. 45, 1-12. 

Suchentrunk, F, Hartl, G. B.; Flux, J. E. C.; 
Parkes,!; Haiden, A.; Tapper, S. (1998): Al¬ 
lozyme heterozygosity and fluctuating asym¬ 
metry in brown hares Lepus europaeus intro¬ 
duced to New Zealand: Developmental 
homeostasis in populations with a bottleneck 
history. Acta Theriol., Suppl. 5, 35-52. 

Swofford D. L.; Selander R. B. (1989): 
BIOSYS-1. A computer program for the ana¬ 
lysis of allelic variation in population genetics 
and biochemical systematics. Release 1.7. 
Users manual. Champaign: Illinois Natural 
History Survey. 

Tiedemann, R.; Hammer, S.; Suchentrunk, F; 
Hartl, G. B. (1997): Allozyme variability in 
medium-sized and large mammals: determi¬ 
nants, estimators, and significance for conser¬ 
vation. Biodiv. Letters 3, 81-91. 

Wright, S. (1978): Evolution and Genetics of Po¬ 
pulations. Vol. IV. Variability within and 
among Natural Populations. Chicago: Univ. 
Chicago Press. 

Yalden, D. W. (1982): When did the mammal fau¬ 
na of the British Isles arrive? Mammal Rev. 
12, 1-56. 

Authors’ addresses: 

Dr. Franz Suchentrunk, Mag. Claudia Jaschke, 

Anita Haiden, Research Institute of Wildlife 

Ecology, Vienna Veterinary University, Savoyen- 

strasse 1, A-1160 Vienna, Austria. 

E-mail: franz.suchentrunk@vv-wien.ac.at 


